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The longest-range part of the nuclear force due to the one-pion exchange governs the energy de¬ 
pendence of the scattering amplitude in the near-threshold region and imposes correlations between 
the coefficients in the effective range expansion. These correlations may be regarded as low-energy 
theorems and are known to hold to a high accuracy in the neutron-proton ^Si partial wave. We 
generalize the low-energy theorems to the case of unphysical pion masses and provide results for the 
correlations between the coefficients in the effective range expansion in this partial wave for pion 
masses up to M-,, ~ 400 MeV. We discuss the implications of our findings for the available and 
upcoming lattice-QCD simulations of two-nucleon observables. 
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I. INTRODUCTION 


The dependence of the nuclear forces and, more generally, nuclear observables upon the fundamental parameters of 
the Standard Model such as the light quark masses is one of the open problems in contemporary theoretical physics. 
Being a fascinating question on its own, it has attracted considerable attention on the theory side and has been 
shown to play important role for various phenomena in nuclear, astro- and particle physics. In particular, Bulgac 
et al. have speculated about an intriguing possibility of the appearance of a bound state in the ^Pq neutron-proton 
partial wave in the so-called chiral limit of Quantum Chromodynamics (QCD) corresponding to massless up and down 
quarks [T]. Such a scenario would clearly have far reaching consequences for nuclear physics. Another interesting 
observation was made in Ref. [2], where the proximity of the physical quark masses to the critical trajectory for an 
infrared renormalization group limit cycle of QCD was pointed out and conjectured to be responsible for the success 
of Efimov’s program for describing the three-nucleon problem, see also Ref. [3] for a follow-up study along this line. 

The quark mass dependence of the nuclear force plays an important role for constraining certain types of extensions 
of the Standard Model. In particular, many such theories allow for the parameters of the Standard Model such as the 
light quark masses to vary over time. The knowledge of the quark mass dependence of the nuclear force combined 
with the theory of Big Bang nucleosynthesis (BBN) and the observed nuclear abundances allow one to constrain a 
possible variation of quark masses at the time of BBN HH]. 

Recently, this topic has experienced renewed interest in connection with the anthropic considerations in nuclear physics 
[5]. In particular, the second O'*' state of ^^C, the so-called Hoyle state, is well-known to play a crucial role for a 
formation of life-essential elements such as ^^C and in hot, old stars due to its closeness to the "‘He-^Be threshold 
[3. Changing its excitation energy of e = 379.47(18) keV, measured with respect to the triple-a threshold, by more 
than ~ 25% was shown to strongly decrease the production rate of either ^^C or BM- The Hoyle state is, 
therefore, viewed as a promising candidate to address the fine-tuning problem of the fundamental constants of Nature 
in connection with the anthropic principle. The sensitivity of the energy difference e to a variation of the light quark 
masses was recently analyzed within an ab initio framework of nuclear lattice simulations mm- This theoretical 
approach makes use of a discretized formulation of chiral effective field theory (EFT) combined with the Monte Carlo 
method to perform Euclidean time evolution of an A-nucleon state, see Ref. m for a review article. Presently, by far 
the dominant source of the theoretical uncertainty in this calculation is related to the lack of knowledge of the quark 
mass dependence of the nuclear force or, more precisely, of the nucleon-nucleon (NN) S-wave scattering lengths. 

While the ultimate answer to the question of quark mass dependence of hadronic observables is eventually to be 
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provided by lattice QCD, one can gain useful insights into this topic within the framework of chiral EFT. During 
the past two decades, this theoretical approach was developed into a powerful tool to derive nuclear two- and many- 
body forces and the corresponding current operators in a systematically improvable and model-independent way, see 
Refs. [131 m] for recent review articles. It exploits the symmetry and symmetry-breaking pattern of QCD to formulate 
an appropriate EFT in terms of pions and nucleons (and possibly of the A(1232) isobar) which are the relevant degrees 
of freedom for low-energy nuclear physics in the non-strange sector. The resulting effective Lagrangian has been used 
to analyze hadronic observables in the Goldstone-boson and single-nucleon sectors as well as nuclear forces and the 
corresponding current operators by means of the chiral expansion, i.e. a perturbative expansion in powers of the pion 
mass Mtt and three-momenta of external particles. In particular, for the two-nucleon force, this expansion has been 
pushed to next-to-next-to-next-to-next-to-leading order [mill! and demonstrated to provide very accurate description 
of low-energy NN scattering observables and the deuteron properties nnin]. Notice that while the number of the 
relevant low-energy constants (LECs) in the effective Lagrangian grows with an increasing order of the calculation, the 
large amount of available experimental data on low-energy proton-proton and neutron-proton scattering observables 
make their determination unproblematic. 

Given its reliance on the spontaneously broken approximate chiral symmetry of QGD, chiral EFT provides, at least 
in principle, a suitable theoretical framework to address quark mass dependence of the nuclear forces in a systematic 
way, see Refs. HHsi nni ns (Mil for recent studies along this line and Refs. [TJ [33] for related calculations at 
the more phenomenological level. Glearly, the major complication which prevents one from being able to perform 
accurate chiral EFT calculations of few- and many-nucleon systems at unphysical values of the quark masses is the 
lack of knowledge of the corresponding LEGs. Quark mass dependence of the short-range interactions induced by 
integrating out the momentum scale associated with real pion production has been discussed in Ref. |26j . Further 
complications emerge due to the appearance of a hnite cutoff in the calculations. This feature is unavoidable if the 
one-pion exchange potential is iterated in the non-relativistic Lippmann-Schwinger equation to all orders. Such an 
approach relies on implicit renormalization which makes it difficult in practice to control the quark mass dependence 
of short-range operators in a systematic way. Notice that the above-mentioned complications can be avoided by 
treating the exchange of pions in perturbation theory using e.g. the approach proposed by Kaplan, Savage and Wise 
(KSW) [27], see Ref. [T9| for studies along this line and Ref. [23] for calculations using a closely related approach with 
dibaryon fields. The perturbative treatment of the one-pion exchange within this framework was, however, shown to 
be inadequate in spin-triplet channels of NN scattering |3H1 [221 • 

Recently, an attempt was made to overcome the above mentioned difficulties associated with the unknown quark mass 
dependence of the short-range part of the nuclear force by employing a resonance saturation hypothesis for contact NN 
operators [3| . It is well known that the LEGs accompanying the NN contact interactions can be understood at a semi- 
quantitative level in terms of the exchange of heavy mesons in the sense of resonance saturation |30j . Using unitarized 
chiral perturbation theory and the information provided by lattice-QGD simulations to determine the pole positions 
of the relevant resonances as functions of the quark masses and assuming the validity of the resonance saturation 
picture at unphysical pion masses allows one to predict the quark mass dependence of the contact operators and to 
carry out chiral extrapolations in the NN sector. This strategy was followed in Ref. |5] up to next-to-next-to-leading 
order in chiral EFT within the Weinberg power counting scheme. It was, in particular, found that the deuteron is 
likely to become less bound at larger values of the light quark mass. These findings are in good agreement with the 
ones achieved in Ref. [23] using a more phenomenological approach by considering a set of Argonne models for NN 
interaction. A detailed comparison of these results with the available earlier chiral EFT calculations both within the 
Weinberg and the KSW power counting schemes can be found in Ref. [3]. We emphasize, however, that modelling 
the Mjr dependence of the contact interaction, as done in Ref. |3] , provides a substantial source of uncertainty which 
was not estimated in that work. 

Finally, chiral extrapolations of the NN S-wave scattering lengths and the deuteron binding energy were considered 
in Ref. [3T] using the chiral EFT formulation of Ref. |32|. In this approach, the leading-order (LO) NN scattering 
amplitude is obtained by solving the three-dimensional integral equation introduced originally by Kadyshevsky [33] . 
This equation is an example of three-dimensional integral equations which satisfy relativistic elastic unitarity. Glearly, 
the Kadyshevsky equation turns into the usual Lippmann-Schwinger equation upon taking the non-relativistic limit of 
the two-nucleon propagator. An important feature of the Kadyshevsky equation for the LO NN scattering amplitude 
is its renormalizability. In particular, all ultraviolet divergences generated by iterations can be explicitly absorbed 
into redefinition of the NN derivative-less contact interaction. This feature allows one to remove the ultraviolet cutoff 
when calculating the amplitude and to avoid the above mentioned complications emerging in the non-relativistic 
framework with a finite cutoff. Assuming that the contributions of the dependent short-range operators are 
suppressed relative to the derivative-less, independent contact terms, the quark mass dependence of the S-wave 
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phase shifts at LO emerges entirely from the propagator in the one-pion exchange potential. The resulting parameter- 
free predictions in the spin-triplet neutron-proton channel are in very good agreement with the calculations reported 
in Ref. [S]. Also the dependence of the ^Sq scattering length agrees qualitatively with the one reported in that 
paper. 

While all above-mentioned studies laiisiisi], which are carried out using rather different theoretical approaches, 
predict less attraction in the deuteron channel for increasing values of the light quark masses, the available lattice-QCD 
simulations seem to indicate an opposite trend. The first pioneering calculation of NN scattering in the framework 
of lattice QCD has been carried out two decades ago by Fukujita et al. [34] within the quenched approximation and 
using pion masses of ~ 550 MeV and heavier. Ten years later, the first fully dynamical lattice QCD calculation of 
the NN S-wave scattering length at pion masses of M,r — 350 MeV, 490 MeV and 590 MeV were reported in Ref. [3S] . 
This study found smaller values of the ^Si scattering length (corresponding to a stronger bound deuteron) but was 
inconclusive of its sign. Notice further that similarly to the calculation reported in Ref. [34], the authors of Ref. [35] 
did not address volume dependence. Recently, binding energies of light nuclei have been calculated in 3-flavor QCD at 
~ 390 MeV [3S] and M^r ~ 800 MeV [37] as well as in 2 -|- 1-flavor QCD at M^r ~ 510 MeV [35], see also Ref. [33] 
for a quenched QCD calculation of the same observables. These fully dynamical lattice QCD calculations found a 
considerably stronger bound deuteron at these large pion masses, see section [TV] for more details. Notice further that 
the NPLQCD Collaboration also extracted the values of the effective range and the first shape coefficient at the pion 
mass of Mt^ ~ 800 MeV [30]. On the other hand, the HAL QCD Collaboration found no bound state in the NN 
^Si-^Di channel by carrying out 3-flavor QCD simulations at pseudoscalar meson masses of 469... 1171 MeV and 
employing a different approach by making use of the two-nucleon potential at the intermediate step of the calculation 
m- Very recently, Yamazaki et al. have performed 2 -|- 1-flavor QCD calculations of light nuclei at pion mass of 
Mtt ~ 300 MeV [33] . They found the deuteron to be even stronger bound at Mt^ = 300 MeV compared to their earlier 
study [39] at — 510 MeV. We further emphasize that the preliminary results of the NPLQCD Collaboration at 
~ 430 MeV [33] also suggest the deuteron to be stronger bound as compared to the physical case. The current 
situation with the lattice-QCD results obtained by different groups is thus not completely clear although there seems 
to be a common trend towards stronger bound deuteron and other light nuclei at values of the light quark masses 
larger than the physical ones. 


In this paper we address the quark mass dependence of low-energy NN scattering observables from the point of view 
of low-energy theorems (LETs). These theorems establish mo del-independent relations between the coefficients in the 
effective range expansion (ERE) which are determined by the long-range part of the nuclear force governed by the 
pion exchange. We generalize the LETs to unphysical pion masses and discuss their accuracy and the validity range 
with respect to the variation in We apply the LETs to provide relations between the coefficients in the ERE in 
the ^Si partial wave for pion masses up to ~ 400 MeV. Our results open the way for nontrivial consistency checks 
of the existing and upcoming lattice-QCD calculations (see also Ref. [33] for a related discussion). As an example of 
a possible application, we employ a linear with respect to extrapolation of the ^Si effective range suggested in 
Ref. [3D] to predict the resulting dependence of the scattering length, shape parameters and the deuteron binding 
energy. The obtained results agree well with the lattice-QCD calculations reported in Refs. [36l|38l|33]. We also use 
the LETs to predict the values of the scattering length, effective range and the shape parameters using the deuteron 
binding energies calculated at = 300 MeV [33], = 390 MeV [3D] and = 430 MeV [33] which can be tested 

on the lattice by calculating the phase shifts. 


Our paper is organized as follows. In section [T^ we discuss in detail the meaning of the LETs and provide results at 
the physical value of the pion mass in the ^Sq and ^Si channels. A g eneralization of the LETs to unphysical values 
of the pion mass in the ^Si partial wave is considered in section III Implications of these findings for lattice-QCD 
results are addressed in section m The main results of our work are summarized in section E] 


II. LOW-ENERGY THEOREMS FOR THE PHYSICAL PION MASS 


Long-range interactions are responsible for the near-threshold left-hand singularities of the partial wave scattering 
amplitude and control its energy dependence [35] [33] . In particular, they impose correlations between the coefficients 
in the effective range expansion which can be regarded as low-energy theorems. In the following, we discuss in some 
detail the meaning of the LETs following the lines of Refs. [35H17] and using the framework of non-relativistic quantum 
mechanics which is appropriate for analyzing low-energy NN scattering. For a related discussion on reconstructing 
the scattering amplitude in the physical region based on the discontinuities across the left-hand cuts and employing 
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FIG. 1: Cut structure of the partial-wave nucleon-nucleon scattering amplitude in the complex energy plane in the absence of 
electromagnetic interactions. The cross refers to a possible bound state. 


the unitarity constraints see Refs. [HI 03] and references therein. 


Consider two non-relativistic particles of mass m interacting via some short-range potential V. The corresponding 
S'-matrix for an uncoupled channel with the orbital angular momentum I is parametrized in terms of a single phase 
shift Si and can be written in terms of the T-matrix as 


Si = = 1- 



( 2 . 1 ) 


where k denotes the scattering momentum in the center-of-mass system (CMS). In the complex energy plane, the 
partial-wave scattering amplitude and thus also the T-matrix possess a so-called unitarity cut, a kinematic singularity 
due to two-body unitarity. The unitarity cut starts from the branch point at the threshold {E = 0) and goes to 
positive infinity. Furthermore, there are singularities associated with the interaction mechanism and located at the 
negative real axis. In particular, in the case of the Yukawa potential (~ exp(—Mr)/r) corresponding to an exchange 
of a meson of mass M, the amplitude has a left-hand cut starting at = —M^/4. Bound and virtual states reside 
as poles at the negative real axis {k = i\k\ and k = —i\k\ for bound- and virtual-state poles, respectively) while 
resonances show up as poles at complex energies. The cut structure of the nucleon-nucleon scattering amplitude in 
the absence of long-range electromagnetic potentials is visualized in Fig. 

It is useful to express the T-matrix in terms of the so-called effective range function T/(fc) = fc^*+^cot(5i(fc) via 


Ti{k) 


TO T;(fc) - ifc2i + l ■ 


( 2 . 2 ) 


Contrary to the scattering amplitude, the effective range function does not possess the kinematic unitarity cut and 
can be shown to be a real meromorphic (i.e. analytic except for poles) function of k^ near the origin k = 0 for 
non-singular potentials of a finite range isniEii- It can, therefore, be Taylor-expanded about the origin leading to the 
well-known effective range expansion (ERE) which has the form 

fc^*''‘^coti5/(fc) = -rfc^ -I- V2k^ + v^k^ + v^k^ -I- ... , (2.3) 

a 2 

where a and r refer to the scattering length and the effective range, respectively, while Vi denote the so-called shape 
parameters. The convergence radius of the ERE is bounded from above by the lowest-lying left-hand singularity 
associated with the potential. In particular, given that the longest-range part of the strong nuclear force is due 
to the one-pion-exchange potential (OPEP), the ERE for NN scattering is expected to converge for energies up to 
l^'iabl ~ M^/(2TOAr) = 10.5 MeV. Notice that the actual convergence range of the ERE might be smaller if the effective 
range function possesses poles corresponding to zeros of the scattering amplitude whose positions are determined by 
the strength of the interaction. Such a situation emerges, for example, when phase shifts change the sign. The problem 
of decreasing the applicability range of the ERE can be easily avoided with if the Taylor expansion is replaced by 
e.g. Pade approximants, see Ref. |52] for a related discussion. Furthermore, we emphasize that a generalization 
of the above considerations to coupled channels such as e.g. the ^Si-^Di channel in neutron-proton scattering is 
straightforward and amounts mainly to replacing the amplitude by a 2 x 2 matrix, see Ref. |53j for more details. 


The framework of the ERE can be generalized to the case in which the potential is given by a sum of long-range 
{tl ~ short-range (rg ~ <C potentials, Vl and Eg, respectively. Following van Haeringen and 
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Kok [^, one can define the modified effective range function via 

1^21 + 1 

F,^{k^) ^ Rf{k) + cot[d,(fc) - 5^{k)]. (2.4) 

In this equation, f["{k) denotes the dost function defined according to f["{k) = ft{^^'^)\r=o ff"{k,r) being the 
dost solution of the Schrodinger equation corresponding to the potential Vl, i.e. the particular solution that fulfils 

lim r) = l. (2.5) 

r—>-oo 

Further, Sl"{k) denotes the phase shift associated with the potential Vl and the quantity Rf(k) can be computed 
from f/"{k,r) as follows: 


R^k) 


lim [ ^i ft{k,r) - 

\ 2 J n r^o [dr2'+i fL{k) 


( 2 . 6 ) 


Here and in what follows, the superscript “L” indicates that the quantity of interest can be calculated solely from 
the long-range part of the potential Vl- As proven in Ref. [S3], the modified effective range function does 

not contain the left-hand singularities associated with the long-range potential and reduces, per construction, to the 
ordinary effective range function Ffik'^) in the case of Vl = 0. It is a real meromorphic function in a much larger 
region set by as compared to Ffik^). In particular, for Yukawa-type potentials, the region in which the modified 
effective range function is meromorphic is set by |fc| < Ms/2. Similarly to the ERE, one can Taylor expand the 
function F^{k‘^) near the origin via 






M, 2 


,.M,4 




fc® + vfk^ 


(2.7) 


This expansion is referred to as the modified effective range expansion (MERE). 


The most frequently used application of the above framework concerns proton-proton scattering. In that case, the 
long-range interaction is due to the Coulomb potential Vl{t) = a/r with a being the fine structure constant. The 
left-hand cut starts directly at the threshold so that the ERE has zero radius of convergence. Notice further that 
the Jost solution and, consequently, the function Rf"{k) can be calculated analytically for the case of the Coulomb 
potential. For example, for I = 0 and the repulsive Coulomb potential, the modified effective range function takes the 
following well-known form: 


F^^ik^) = Fc{k^) = C^ir,) k cot[Sik) - S^ik)] + 2kr]hir ]), 


( 2 . 8 ) 


where the Coulomb phase shift is = arg r(l + irf) and the quantity rj is given by 


^=Tk^- 

Further, the functions (^ 0 ( 77 ) (the Sommerfeld factor) and h{r]) read 


- g27rr; _ ^ > 


and h{r]) = Re 47( 777 ) — In)?;). 


(2.9) 


( 2 . 10 ) 


Here, 4'(z) = V{z)/V{z) denotes the digamma function. For more details on the analytic properties of the scattering 
amplitude and related topics the reader is referred to the review article [55j . 

We are now in the position to clarify the meaning of the LETs. Assuming that the functions Ffikfi) and F^{kfi) do 
not possess discrete poles in their meromorphic regions, the size of the coefficients in the ERE and MERE (except 
for the scattering length) is expected to be governed by the scales Ml and Ms associated with the lowest left-hand 
singularities, see |33| for a related discussio n. I f the long-range interaction is known, the quantities //"(fc), Rf{k) and 
i5/^(A:) entering the right-hand side of Eq. (|2.4|) can be calculat ed ex plicitly. Using the MERE, truncated at a given 


order to approximate the function F^{k‘^)^ and inverting Eq. (|2.4|) then allows one to express the phase shift 5i{k) 
in terms of the known quantities fi"{k), Rf{k) and 5l{k) and the first coefficients in the MERE which parametrize 
physics associated with the short-range interaction. For example, at LO, the modified effective range function is 
approximated as F^(k^) ~ Thus, using a single piece of information about the short-range interaction in 
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form of or, equivalently, the usual scattering length^ a allows one to predict all coefficients in the ERE. These 
predictions are accurate up to corrections emerging from the second term in the MERE. The MERE for {k'^) 
thus provides a systematically improvable expansion of the subthreshold parameters entering Eq. ( |2.3[ ) in powers of 
Ml/M s- The resulting correlations between the subthreshold parameters at each fixed order in this expansion are 
referred to as LETs. For an example of an exactly solvable model, which allows one to get analytical insights into the 
meaning of the LETs, the reader is referred to Ref. [45] . 

We are now in the position to address the LETs for S-wave neutron-proton scattering. Here and in what follows, 
we restrict ourselves to the case of exact isospin symmetry and neglect the long-range electromagnetic interactions 
between the neutron and proton. The longest-range part of the nuclear force is due to the OPEP which in momentum 
space has the form 


Vi.(g) 


g\ ai-qa 2 -q 
4E2 q^+M^ ^1-^2, 


( 2 . 11 ) 


where q = p' — p is the momentum transfer of the nucleon while cfj (r^) denote the spin (isospin) Pauli matrices of the 
nucleon i. Further, gA and EV refer to the axial vector coupling of the nucleon and pion decay constant, respectively. 
Throughout this work, we adopt the values of = 92.4 MeV, gA = 1.267 for these constants at the physical value 
of the pion mass = (2M,r± + -^7ro)/3 = 138.03 MeV. Upon performing Fourier transform, one obtains the 

well-known coordinate-space expression for the long-range part of the OPEP, namely 


Ui^(r) 


M A 

127r 4E2 


T 2 






+ (Tl ■ <72 


^-M„r 

M^r 


( 2 . 12 ) 


where S 12 is the usual tensor operator, S 12 = 3(cti • 'r){<72 ■ r )/r^ — ai- $ 2 - Notice that the above potential describing 
point-like nucleons interacting via one-pion exchange is singular and thus cannot be directly employed in the calcu¬ 
lation of the LETs using the quantum mechanical framework outlined above.^ One possible way to overcome this 
complication is to apply a properly regularized version of the OPEP. For example, one can employ a coordinate-space 
regularization procedure introduced in Refs. ISaEH] or a closely related method used in Refs. imiiTi. These ap¬ 
proaches are analogous to the Pauli-Villars regularization and maintain, per construction, the analytic structure of 
the amplitude in the sense that non-analyticities induced by the regulator appear at momenta k > A. Here, A refers 
to the momentum scale associated with the regulator which in such an approach should be chosen of the order of 
the breakdown scale in the problem 03 Eg. The properly regularized OPEP can then be used to compute the dost 
solution by solving the Schrodinger equation and to work out the LETs for NN scattering. 


Instead of using the approach described above, we employ here the modified version of Weinberg’s chiral EFT proposed 
in Ref. [32] in order to work out the LETs for NN scattering. This momentum-space framework is more convenient 
for performing numerical calculations of the LETs and, as already mentioned in the introduction, has an advantage 
of being well suited for carrying out chiral extrapolations beyond the LETs |5D] . 

Within the modified Weinberg EFT approach of Ref. 02], the LO amplitude is obtained by solving the Kadyshevsky 
equation 03] which, for the case of the fully off-shell kinematics, takes the form 


T ipo,p\p) = V {p\p) + 


(fqV{p',q) G{po,q) T{po,q,p) , 


where G{pQ,q) is the free Green function 


G{pQ,q) 


m% _ 1 _ 

2(2 7r)3 ^-2 _ y^q2 ^ 


(2.13) 


(2.14) 


^ Note that is directly related to the scattering length a via Eq. | |2.4[ |. Similarly, the modified effective range can be calculated 

using a and r as input. _ _ 

^ It should, however, be understood that the expressions for the OPEP in Eqs. ( |2.1l[ |, ( |2.12[ | become meaningless at short distances 
where effects due to finite nucleon size cannot be neglected, see Ref. ESI for a related discussion. Furthermore, while the OPEP clearly 
dominates the interaction between the nucleons at large distances, the two-pion exchange potential becomes comparable in size to the 
OPEP at distances smaller than r ~ 2 fm. 
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FIG. 2: Correlations between the inverse scattering length a“^, effective range r and the first three shape parameters V 2 , V 3 
and V 4 induced by the one-pion exchange interaction in the ^So (left panel) and ^Si (right panel) channels. Solid rectangles, 
diamonds, open triangles and circles correspond to the values of r, V 2 , vs and V 4 , respectively, extracted from the Nijmegen 
partial wave analysis [Sam]. The results in the ®Si partial wave correspond to the Blatt-Biedenharn parametrization of the 
S-matrix |62| . 


Further, p {p') is the incoming (outgoing) three-momentum of the nucleon in the center-of-mass frame, po = 

with iriN denoting the nucleon mass and k being the corresponding three-momentum of an incoming 
(on-mass-shell) nucleon. The integral equation (2.13) is solved numerically in the partial wave basis. 


In this work, we are interested in the ^Sq and the coupled ^Si-^Di channels of neutron-proton scattering. In both 
cases, the LO potential is given by the OPEP specified in Eq. (2.11) accompanied by a momentum-independent 


contact interaction, whose strength is adjusted to reproduce the S-wave scattering length. We regularize the integral 
equation by introducing an ultraviolet momentum cutoff A and take the limit A —>■ 00 when calculating the phase 
shifts. Notice that contrary to the non-relativistic framework based on the Lippmann-Schwinger equation, this is a 
legitimate procedure since the LO integral equation is renormalizable in the sense that all divergences emerging from 
its iterations are absorbable into redefinition of the LO contact interaction. Once the phase shifts are calculated, we 
extract the effective range and the shape parameters numerically by matching the truncated ERE to the effective range 
function. Given that the calculated S-wave scattering amplitudes fulfill elastic unitarity, correctly reproduce the first 
left-hand cut due to the OPEP and match the experimental value at the threshold (given by the scattering length), 
the resulting predictions for r and Vi are equivalent to the LO LETs discussed above in the context of the MERE. 
We emphasize, however, that employing the cutoff-independent framework of Ref. |32j is not crucial for testing the 
LETs in NN scattering. A more conventional non-relativistic chiral EFT approach based on the Lippmann-Schwinger 
equation and utilizing a finite ultraviolet cutoff is equally well suited for this purpose. 


The LO LETs for the ^Sq and the ^Si channels have already been addressed in Ref. [35]. In Eig. we show the 
LO LET predictions for the effective range and the first three shape parameters as functions of the inverse scattering 
length which is an input parameter in our calculations. Here and in what follows, the results in the ^Si partial wave 
correspond to the Blatt-Biedenharn parametrization of the S-matrix |62] . To avoid the appearance of large/small 
numbers, we use dimensionless quantities by multiplying the coefficients in the ERE with the corresponding powers 
of the pion mass. As expected, the predicted dimensionless coefficients appear to be of a natural size for a rather 
large range of values of the scattering length. Notice further that the limit a —>■ ± 00 , which describes the situation 
in which there is a zero-energy bound state, does not lead to any peculiarities in r and Vi. In both channels, one 
observes a clear tendency for r and Vi to become unnaturally large in magnitude in the regions of l/(aM,r) ^ —0.5 and 
l/(aM 7 r) > 1. We found that the unnaturally large values of r and Vi for l/(aM.n.) > 1 are due to the appearance of 
a pole in the effective range function at relatively low positive energy. The pole position corresponds to the energy at 
which the corresponding phase shift crosses zero. Eor example, for l/(aM,n.) ~ 1.4, the pole in the ^Si effective range 










TABLE I: Various available results for the low-energy theorems for the neutron-proton ^So and ^Si partial waves as described 
in the text. The results in the ^Si partial wave correspond to the Blatt-Biedenharn parametrization of the S-matrix | 62| . 



a [fm] 

r [fm] 

V 2 [fm®] 

V 3 [fm®] 

V 4 [fm^j 

Neutron-proton ^Sq partial wave 
LO, Ref. [31] 

fit 

1.50 

-1.9 

8.6(8) 

-37(10) 

NLO, nonperturbative C 2 , Ref. [63] 

fit 

fit 

-0.61... -0.55 

5.1... 5.5 

-30.8... - 29.6 

NLO KSW, Ref. [H] 

fit 

fit 

-3.3 

18 

-108 

Empirical values, Ref. [53] 

-23.7 

2.67 

-0.5 

4.0 

-20 

Nentron-proton ®Si partial wave 
LO, Ref. [31] 

fit 

1.60 

-0.05 

0.82 

-5.0 

NLO, this work 

fit 

fit 

0.06 

0.70 

-4.0 

NLO KSW, Ref. [H] 

fit 

fit 

-0.95 

4.6 

-25 

Empirical values, Ref. [64] 

5.42 

1.75 

0.04 

0.67 

-4.0 


function is located at A: ~ 60 MeV which is already within its meromorphic region. For larger values of l/(aM 7 r), the 
pole moves towards the threshold so that the values of the r and vi are actually governed by the pole position rather 
than by M^. As already emphasized above, the appearance of a pole in the effective range function can be easily 
accommodated by replacing the Taylor expansion of by e.g. Fade approximation. On the other hand, the case 

l/(aM 7 r) < —0.5 corresponds to the absence of a bound state and, given the natural value of the scattering length, 
it might already describe a perturbative regime of NN scattering. For weakly interacting systems, it is easy to see 
analytically (by making use of the Born approximation) that the coefficients in the ERE do not scale anymore with 
the range of the interaction. This is due to the appearance of an additional dimensionless small parameter associated 
with the weakness of the interaction. 


The predictions for r and Vi based on the LO LETs and corresponding to the experimentally observed values of the 
scattering lengths are compared in Fig. with the empirical numbers from Ref. |531l61j based on the Nijmegen partial 
wave analysis (PWA). We also collect various available results for the LETs in the ^Sq and the ^Si channels in Table 
|Tj The results for the LO LETs shown in the table are taken from Ref. [32]. While the predicted values of r and Vi 
do reproduce the empirical values in the ^Si channel rather accurately, the agreement in the ^Sq partial wave is, at 
best, on a qualitative level. This pattern can be understood as follows. Generally, the accuracy of the LO LET is 
expected to be set by the ratio El/Es, where El (Es) denotes the energy corresponding to the branch point of the 
left-hand cut due to the long-range interaction (branch point of the first left-hand cut in the scattering amplitude 
which is not correctly described by the employed approximation for the interaction). Therefore, given that we do not 
take into account the two-pion exchange potential in our calculations and thus cannot describe correctly the second 
left-hand cut, one may expect the LO LETs to be accurate at the level of ~ 25%. In fact, the accuracy of the LO 
LETs in the ^Si partial wave appears even to be somewhat higher.^ This is presumably due to the fact, that we 
actually include the parts of the second and higher left hand cuts (see Fig. which are associated with iterations of 
the OPEP. In fact, according to the chiral power counting, such iterative contributions to the scattering amplitude are 
expected to be more important than the irreducible two- and more-pion exchange potentials. On the other hand, the 
low accuracy of the LETs in the ^Sq channel is a consequence of the weakness of the OPEP projected on that partial 
wave. The strong tensor part of the OPEP, which accompanies the operator S '12 in Eq. (2.12), does not contribute 
to spin-singlet partial waves. The weakness of the OPEP in the ^Sq channel is evidenced by the fact, that the phase 
shift corresponding to the long-range OPEP reaches only about 12.6° at maximum, which has to be compared with 
max{Siso(k)) ^ 65° in the real world. 


We also list in Table |Tj for the sake of completeness, the predictions obtained within the KSW framework [27] 
which relies upon a perturbative treatment of the OPEP. These results are taken from Ref. [55] and show very large 
deviations from the empirical values. This shows that the perturbative treatment of the OPEP in these channels is 


® The large relative deviation for the first shape parameter should not be taken too seriously due to its unnaturally small value. 
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not appropriate. 

In Ref. [53], the subleading contact interaction C 2 {p^ +p'‘^) has been taken into account in the ^Sq channel within the 
framework of Ref. |32] , both perturbatively and nonperturbatively. Given that the long-range interaction included in 
that work is entirely given by the OPEP, the resulting predictions for the shape parameters, which are listed in Table 
[I] for the case of the nonperturbative treatment of C 2 , are equivalent to the next-to-leading-order (NLO) LETs. Notice 
that the non-perturbative inclusion of the subleading contact interaction necessarily results in a residual subtraction 
scale dependence of the amplitude which is why the predictions for ^ 2 , 3,4 are shown within a range of values, see 
Ref. [63] for more details. While a significant improvement is observed in the description of all Vi, the convergence 
is clearly rather slow in this channel. We find that the accuracy of the LETs in the ^Sq channel is insufficient for 
the purpose of providing constraints on the available lattice QCD results. Eor this reason, we will concentrate in the 
following entirely on the spin-triplet channel. 

While the LETs in the ^Si channel appear to be fairly accurate already at LO as far as one is concerned with predicting 
the values of r and Vi from the scattering length, their accuracy is still insufficient for extracting the deuteron binding 
energy using the effective range as input, as we intend to do in the next sections. This is because of the function r{a~^) 
being rather flat in the relevant region of values of l/(aM^), which tends to magnify the uncertainty in the value of a~^ 
extracted from r. Notice that the observed nearly quadratic dependence of r on a~^ may be related to the universal 
behavior of this function found for van der Waals-like interactions |65|, see also Ref. ^6] for a related discussion. In 
addition, the deuteron binding energy depends, to a good accuracy, quadratically on a~^, so that the error in the 
value of the binding energy is further magnified. For example, using the LO LETs together with the experimentally 
observed value of the effective range r = 1.75 fm as input, one extracts a = 7.2 fm and B 4 = l.I MeV which is two 
times smaller than the experimentally observed value of Bd — 2.22 MeV. To deal with this issue, we extend the LETs 
in the ^Si channel to NLO by taking into account a subleading short-range interaction whose strength is adjusted to 
reproduce the experimentally observed value of the effective range. Notice that contrary to the ^Sq channel, where the 
long-range part of the OPEP is non-singular and a non-perturbative inclusion of the subleading contact interaction can 
be carried out semi-analytically in a close analogy to pionless EFT, see Ref. [63| for more details, we are unable to treat 
the subleading contact interactions non-perturbatively in spin-triplet channels without, at the same time, destroying 
the explicit renormalizability feature of the scattering amplitude. A perturbative inclusion of the NLO short-range 
terms will be discussed in detail in a subsequent publication EOj. Here, in order to avoid possible shortcomings due to 
the perturbative treatment of the subleading short-range terms, we follow a different approach and employ resonance 
saturation to model higher-order contact interactions by means of a heavy-meson-exchange. In particular, we use for 
the NLO potential 


Rnlo = kill-( 9 ) + Co + (3 


ai- qa 2 - q 


(2.15) 


where the heavy-meson mass M is set M = 700 MeV and the strength f3 is adjusted to reproduce the empirical value 
of the effective range. In agreement with the arguments given above, we have verified that the results for the LETs 
are insensitive to details of the short-range interaction. In particular, the deviation in the results caused by using 
the sigma-like scalar potential instead of the tensor one in Eq. (2.15) for the sub-leading short-range interaction is 
negligibly small. The sub-leading short-range interaction is only needed to describe the second term in the MERE 
^ which appears, in addition to a^, as an input parameter in our calculation. The employed form of the NLO 
potential ensures that the corresponding scattering amplitude is renormalizable so that the NLO calculations can be 
carried out in the same way as at LO. The predicted values of the shape parameters after fixing the values of Co and 
/3 are listed in Table [l| and, as expected, show a clear improvement as compared with the LO results. In fact, the 
NLO LETs appear to be accurate at the level of a few percent (except for V 2 ). 


III. LOW-ENERGY THEOREMS FOR UNPHYSICAL PION MASSES 


The LETs discussed in the previous section can be straightforwardly generalized to the case of non-physical quark 
(or pion) masses. The main effect from changing the pion mass, as far as the analytic structure of the amplitude is 
concerned, corresponds to shifts of the branch points of the left-hand cuts due to one-, two- and more-pion exchange, 
see Fig. [2 This effect originates from the explicit pion mass dependence in the pion propagator. Next, one needs to 
account for the change in the discontinuity across the left-hand cuts caused by the dependence of the strength 
of the OPEP, namely the ratio qaIFtt. To account for this effect, we make use of the lattice-QCD results for the 
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FIG. 3: Quadratic polynomial regression fits to lattice QCD data for the pion decay constant fV, nucleon mass niM and the 
nucleon axial-vector coupling constant qa- Lattice data for FV and rriM (qa) correspond to the 2-1-1 flavor simulations of 
the Budapest-Marseille-Wuppertal collaboration using Wilson fermions m and 2 and 2 -|- 1 -|- 1 flavor simulations reported 
in Ref. |68] (Ref. [^) using twisted mass fermion ensembles. For details of the simulations and analyses see the above 
mentioned references. Filled circles without error bars show the experimentally measured values. The dashed lines depict our 
fit results to the lattice data and the experimental values while the shaded bands correspond to the 67% confidence levels of 
the interpolations. 


dependence of these quantities, see Fig. In particular, we performed quadratic polynomial regression fits 
(as functions of Mj) of the results from the 2-1-1 flavor simulations by the Budapest-Marseille-Wuppertal (BMW) 
collaboration for m and from the 2 and 2 -|- 1 -|- 1 flavor simulations reported in Ref. [55] for qa, see also Ref. [70] 
for recent lattice-QCD calculations of this quantity. In our fits, we included lattice-QCD data for pion masses up to 
= 500 MeV. Further, we enforced the fits to go through the experimental value for a given quantity at the physical 
value of the pion mass. The results of our fits together with the 67% confidence region are shown in Fig. Notice 
that the observed tendency of lattice-QCD results for the axial charge of the nucleon to somewhat underestimate the 
experimental value may indicate that certain systematic corrections have not been properly taken into account in 
these studies. We expect this issue to be clarified in the near future, when more results near the physical point will 
become available and the statistical uncertainties of the simulations will be improved. Notice further that we did not 
include all available lattice QCD results in our fits since this will unlikely improve the accuracy of the interpolations. 
Especially for the pion decay constant, the uncertainty of the interpolation of the lattice-QCD results is very small 
and can be safely neglected at the accuracy level of our calculations. 

Next, one also has to account for the pion mass dependence of the nucleon mass. Here, we follow the same strategy 
as for Ft^ and qa and employ a quadratic interpolation of the recent lattice-QCD results for mN{M^) of the BMW 
collaboration [57], see Fig. We emphasize that our fits of the dependence of gA and are not intended 
to provide correct chiral extrapolations of these quantities to the chiral limit but serve exclusively to enable a smooth 
interpolation of the available lattice-QCD results. 

Using the above results for the pion mass dependence of F^., gA and uin together with the explicit M^. dependence 
in the pion propagator allows us to describe the discontinuity across the first left-hand cuts due to the OPEP for 
arbitrary values of the pion masses. Thus, we can immediately generalize the LO LETs by calculating r and Vi as 
functions of the (inverse) scattering length at unphysical values of the pion mass. The resulting predictions are shown 
by the various thick lines in Fig. [^ for = 50, 100, 150, 200, 300 and 400 MeV. Notice that contrary to chiral EFT 
extrapolations, we do not make here any assumptions about the quark mass dependence of the included derivative-less 
short-range interaction cx Cq, which is traded into the dependence of the calculated observables on the value of the 
scattering length. In particular, we perform calculations with different pion masses as if we lived in different worlds, 
in which the physical value of the pion mass were M,r = 50, 100, ... MeV. Thus for each given pion mass the value 
of Co is adjusted to reproduce the given value of the scattering length used as input. In this sense, our results are 
more general than the ones obtained for the ^Sq and ^Si channels in Ref. m at LO within the modified Weinberg 
approach, where the dependence of the short-range interaction on was neglected in agreement with the estimation 
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l/{aM„) 




l/{aM^) 





FIG. 4: Correlations between the inverse scattering length , effective range r and the first three shape parameters V 2 , V 3 
and V 4 in the ®Si partial wave induced by the one-pion exchange interaction. Various thick lines show the predictions of the 
LO LETs while light-shaded bands (hardly visible for small M-^) between thin lines depict the results of NLO LETs and reflect 
the estimated uncertainty due to unknown M-^ dependence of the subleading short-range interaction as explained in the text. 


based on naive dimensional analysis (NDA). 


In order to extend these results to NLO, we need to specify the M ^r dependence of the subleading short-range 
interaction which is modeled via resonance saturation, see Eq. (2.151. While this effect would be suppressed in the 


chiral EFT approach if one assumes the validity of NDA for short-range interactions, we decided to estimate it in 
order to stay as model independent as possible. Specifically, we assume that the dependence of the strength f3 of 
the short-range interaction is within the envelope built by the lines which go through the physical point and describe 
a ±50% change in the value of /3 for = 500 MeV, i.e.: 


i-sp 


Ml - 


AM2 


< 


P{M^) 


<i + sp 


Ml - 


AM2 


(3.16) 


with SP = 0.5 and AMl = {Ml — (AfP*'y®)2)|^ MeV' estimation is justified for the considered 

quantities such as gA, Ett and rriN, as can be seen from Fig. 

Our NLO LETs predictions for the effective range and the shape parameters viewed as functions of the inverse 
scattering length are visualized in Fig. [4] by th e light-shaded bands. These bands correspond to the variation of P 
at a given value of M^. according to Eq. ( |3.16 ). For pion masses below M^^ ~ 200 MeV, one observes a rather small 
difference between the predictions based on LO and NLO LETs whose size can be viewed as an estimation of the 
accuracy of the LO LETs. Also, the employed variation of P is essentially invisible for such values of the pion mass. 
For heavier pions, both the differences between the LO and NLO LETs as well as the uncertainty in the results 
associated with the Mt^ dependence of P start to increase. As shown in this figure, one cannot expect to have accurate 
predictions for pion masses above M^^ ~ 400 MeV. Notice that decrease in the predictive power of the LETs is to be 
expected for heavier pions due to the decreasing separation between the soft and hard scales in the problem. 
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TABLE II: Available experimental and infinite-volume lattice QCD data for nucleon-nucleon scattering parameters and bound 
state energies in the ^So and ^Si channels at various values of the pion mass. 



= 138 MeV 

= 300 MeV [12] 

Mtt = 390 MeV [35] 

Mtt = 510 MeV [33] 

= 800 MeV [ID] 

The ®Si channel 
Bd [MeV] 
a [fm] 
r [fm] 

2.224 

5.42 

1.75 

14.5(0.7)(1^;^) 
not given 

not given 

11(05)(12) 

not given 

not given 

11.5(1.1)(0.6) 

not given 

not given 

19.5(3.6)(3.1)(0.2) 

1 . 82 («; 1 ^)(«: 1 D 

0.906(irof5)(io.Sf4) 

The ^So channel 
[MeV] 
a [fm] 
r [fm] 

-23.7 

2.67 

8.5(0.7)(l^j) 
not given 

not given 

7.1(5.2)(7.3) 

not given 

not given 

7.4(1.3)(0.6) 

not given 

not given 

15.9(2.7)(2.7)(0.2) 

9 qQ/+0.19\/+0.27\ 

1 iQn/+0.071w+0.059\ 
077/1-0.063 J 




FIG. 5: Left panel: Linear with interpolation of the quantity M-^r in the ^Si partial wave according to Eq. (4.171 as 
sugg ested in Ref. m- Right panel: Linear with interpolation of the quantity M^r in the ®Si partial wave according to 
Eq. (4.181. In both cases, solid squares refer to the value of at the physical point. 


IV. IMPLICATIONS FOR LATTICE QCD CALCULATIONS 


We are now in the position to confront the LETs with the available lattice QCD results in the NN sector. In table 
|TTj we list the published lattice-QCD results for the S-wave scattering parameters and energies of the bound states 
together with the experimental data. We do not show in the table the results from Ref. [5S] where volume dependence 
was not addressed. Unfortunately, lattice calculations in the NN sector focused so far mainly on the binding energies 
and do not provide information on the scattering parameters. An exception is the work of Ref. [40) . which provides, 
in addition to the binding energies, also the values of the scattering length, effective range and even the first shape 
parameter at the pion mass of M,r — 800 MeV. Clearly, such heavy pion masses are out of reach of the LETs. On the 
other hand, the authors of Ref. [40) conjectured that the quantity may exhibit a nearly linear dependence on 
the pion mass. The suggested linear interpolation between the physical point and the lattice result has the form [40) : 


M^r - where = 0.7261^^5^1^59 , 


= 3.701°;^^!°;^^ GeV^, 


(4.17) 


and is visualized in the left panel of Fig. While we cannot judge on the validity of the suggested linear dependence 
of the quantity on the pion mass based on the LETs alone, we can test its compatibility with the lattice-QCD 
results for the deuteron binding energy available for pion masses within the validity range of the LETs. Specifically, 
we employ the effective range r(Mn.) from Eq. (4.17) instead of the scattering length to fix the dependence of 


the short-range interaction by adjusting the value of Cq in Eq. (2.15) and make predictions for the scattering length. 
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FIG. 6 : NLO LET predictions for the pion mass dependence of the deuteron binding energy, the ratio 'yd/M-n, the ratio a/r 
and the first three shape parameters in the ®Si partial wave assuming the linear M„ dependence of the effective range specihed 
in Eq. (4.171 and visualized in the left panel of Fig. ^ Dark-shaded bands show our estimation of the uncertainty of the NLO 
LETs due to the unknown M-^ dependence of the subleading short-range interaction specified in Eq. (3.161, light-shaded bands 
depict the uncertainty in the linear extrapolation of the effective range used as input, as shown in the left panel of Fig. 


shape parameters and the deuteron binding energy. Our results for the deuteron binding energy Bd, the ratio 
where 7 ^ = y/Bdmyy is the deuteron binding momentum, the ratio a/r and the first three shape parameters M^V 2 , 
M^vs and MIv 4 are visualized in Fig. In this figure, the dark shaded bands result from the variation of the 
constant /3 specified in Eq. (3.16) and reflect the uncertainty of the NLO LETs^. The light-shaded bands correspond 


to the resulting uncertainty which emerges from the theoretical uncertainty at NLO and the errors of the linear 
interpolation of M^r{MT^) (see the left panel of Fig. [5 and Eq. ( 4.17| )) added in quadrature. Notice that we also show 
in Fig. 1^ the preliminary lattice-QCD result of the NPLQCD collaboration at M,r = 430 MeV |33]. Remarkably, the 
linear dependence of suggested in Ref. m indeed appears to describe very well the common trend of the 
lattice-QCD results for the deuteron binding energy at intermediate pion masses. Also the NPLQCD Collaboration 
results of Ref. m for Bd, 7 d/M^, a/r and M^V 2 at the pion mass of = 800 MeV can be well described by further 
extrapolating our results to heavier pion masses without introducing any strong curvature. Assuming the validity of 
Eq. (4.17) for pion masses below the physical one, we conclude that the deuteron becomes unbound for M^r ~ 50 MeV. 


It is also interesting to notice that the scattering length and the shape parameters show rather strong variations with 
the pion mass around and below the physical point. This nontrivial behavior is driven by the long-range physics 
associated with the pion exchange and is, in principle, testable in lattice QCD. The obtained results for the quantities 
7 d/M.n. and a/r, which probe the amount of fine tuning in the NN system, suggest that the physically realized value 
of the quark mass is close to the point, which separates the strong fine-tuning regime characterized by the rapidly 
growing scattering length from the regime featuring a fairly small amount of fine tuning with a/r = 2... 3 within the 
large range of pion masses. 


We emphasize that the observed agreement between the predicted dependence of the deuteron binding energy and 


^ Note that employing the scalar sub-leading potential in Eq. ( |2.15| l instead of the tensor one yields the results which are well within the 
dark shaded band for all quantities except for the parameter which is relatively small and appears to be slightly outside of this band 
for > 250 MeV. 
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FIG. 7: NLO LET predictions for the pion mass dependence of the deuteron binding energy, the ratio the ratio a/r 

and the first three shape parameters in the ®Si partial wave assuming the linear dependence of the effective range specihed 
in Eq. (4.181 and visualized in the right panel of Fig. For notations see Fig. 


lattice-QCD results is a rather nontrivial consequence of the assumed extrapolation of M^r. To illustrate this point, 
we consider an alternative scenario by assuming that M^r is a linear function of rather than of M^. Notice that 
while a linear in dependence seems to be natural for NN observables, such a scenario cannot be excluded a priori. 
Using the lattice-QCD result of Ref. [ID] for M^r at the pion mass of = 800 MeV, the resulting interpolation 
formula takes the form 


M^r ^ -b where 


ri^Si) _ -I - 140 + 0 . 009 + 0.011 
~ -‘-•'‘-^^-0.009-0.009 ’ 


nC^'S'i) _ q Qr+O.45+0.45 

u — 0.yj_Q_4g_Q_55 


GeV 


-2 


(4.18) 


see the right panel of Fig. The resulting LET predictions for the deuteron binding energy are shown in Fig.[^ The 
assumed quadratic dependence on is clearly in conflict with the lattice-QCD results and, therefore, appears to 
be highly unlikely. Interestingly, the shape parameters seem to be more robust with respect to the variation of the 
functional form of r(M^) and show a qualitatively similar behavior in both scenarios. 


Finally, as a last application of the LETs, we use the quoted values of the deuteron binding energy to predict the 
scattering length, effective range and the first three shape parameters at pion masses of M,r = 300 MeV, M,r = 
390 MeV and M.^. = 430 MeV where lattice-QCD results are already available. Our results for the coefficients in the 
ERE are collected in Table |IlI| In the case of the preliminary NPLQCD calculation at = 430 MeV, the values for 
the deuteron binding energy and the errors are extracted from a corresponding plot of Ref. |43j . These results may 
serve as useful consistency checks for the already published lattice-QCD results if the corresponding lattice data can 
be used to extract, in addition to the deuteron binding energy, also information about the scattering parameters. 


V. SUMMARY AND CONCLUSIONS 


The quark mass dependence of the low-energy NN scattering observables is investigated using the low-energy theorems 
which establish model independent relations between the coefficients in the effective range expansion governed by the 
long-range one-pion exchange potential. In order to clarify the meaning of the LETs, we employed the modified 
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TABLE III: Predictions for the scattering length, effective range and shape parameters based on the lattice-QCD results for 
the deuteron binding energy. The first two errors correspond to the ones of the lattice-QCD results for the binding energy. The 
last quoted error is our estimation of the uncertainty of the NLO LETs associated with the unknown dependence of the 
subleading short-range interaction as described in the text. 



= 300 MeV [42] 

= 390 MeV [36] 

= 430 MeV [33] 

aM-n 

3.59(^°°D(1°o:?D(0.03) 

4 QO/+1.21W+00 W+0.07N 

60/1-1.09/v-0.09/ 

A C9/+0.19w+0.36\/+0.10\ 
"±-'9Z[^_Qj^y;[^_0_27A-0.14J 

tMt, 

1.70(0.01)(ini)(0.03) 

9 99/+0.20w+1.13\/+0.08n 
Z.ZZ[^_0.i2A-0.25/\-0.10/ 

9 97/+0.04w+0.08\/+0.10\ 

' l-O.O 5 A-O.O 7 A-O.i 4 J 

V2MI 

0.23(0.01)(0.01)(1H?) 

r. r,^/+0.06w+0.24w+0.07x 
U.OU[^_0.05A-0.10/\-0.09/ 

A r 0 /+0.02 \/-I-0.03 \ /+0.10 \ 
U.OC>[^_0.o2JI-0.03J\-0.13J 

V3MI 

0 .11(0.01)(0.01)(0.01) 

'J-^'JV_0.09A-0.16A-0.07/ 

0.46(«°^)(«:“)(0.12) 

ViMl 

-0.04(0.01)(1H?)(0-01) 

A riO(+0.05\f+0.09\f+0.11\ 

'-'•'-'OV_O. 13 A-O. 33 A-O.IO/ 

rj -|/i/+0.00\/'+0.00w+0.15\ 

'J--L9:[^_0.oi J 1-0.01/v-0.16/ 


effective range expansion [54] to parametrize the short-range physics in a systematic way. Using the knowledge of the 
long-range interaction together with the MERE truncated at a given order allows one to calculate all coefficients of 
the ERE within a systematically improvable expansion in powers of the ratio of the long- and short-range scales. The 
explicit treatment of the left-hand cuts in the partial-wave scattering amplitudes due to the long-range pionic forces 
results in correlations between the coefficients in the ERE which are regarded as the LETs. 

In this work we considered the LETs for NN scattering in the and channels. At leading order in the above 
mentioned expansion, the short-range physics is parameterized solely in terms of the modified scattering length. 
Therefore, in order to obtain LETs at LO, we considered two nucleons interacting via the OPEP supplemented by 
two contact interactions without derivatives, whose strengths were adjusted to reproduce the scattering length in 
each of these channels. The NN scattering amplitude was calculated by solving the Kadyshevsky equation which is 
exactly renormalizable at this order [52] ■ In the ^Si channel, the predicted values of the effective range and the shape 
parameters agree rather well with the empirical values extracted from the Nijmegen PWA [Mj already at LO. On the 
other hand, the OPEP projected onto the ^5'o partial wave is very weak, so that the LETs show much less predictive 
power in that channel. We, therefore, restricted ourselves to the spin-triplet channel in this paper. To further increase 
the accuracy of our predictions, we extended the calculations to NLO LETs by including the subleading short-range 
interaction and tuning its strength to reproduce the effective range. The resulting predictions for the shape parameters 
in the partial wave are found to be accurate at the level of a few percent. 

As a next step, we generalized the LETs in the ^Si channel to study correlations between various low-energy observ¬ 
ables at unphysical pion masses up to Mt^ — 400 MeV. To this aim, we made use of the available lattice-QCD results 
for the pion mass dependence of the pion decay constant, nucleon axial vector coupling and the nucleon mass. The 
predicted correlations between the low-energy NN parameters at unphysical pion masses open the way for nontrivial 
consistency checks of the ongoing and upcoming lattice QCD calculations in the NN sector and might also be useful 
for reducing the systematic uncertainty. 

As an application, we extracted the low-energy parameters, namely the scattering length, effective range and first three 
shape parameters, from the lattice-QCD calculations of the deuteron binding energy at the pion masses M,r = 300 
MeV [32] j 390 MeV |5^ and 430 MeV |35j. Our predictions illustrate the importance of a simultaneous extraction of 
the scattering parameters in addition to the binding energies in lattice-QCD simulations of the NN system. 

As another application, we applied the LETs to predict the behavior of the binding energy, scattering length and 
shape parameters based on the linear dependence of the effective range suggested in Ref. |30|. We found 
that the predicted shape of the binding energy as a function of M.^. is in good agreement with the general trend of the 
lattice QCD data. We also found that the scattering length and the shape parameters show rather strong variations 
with the pion mass around and below the physical point. This nontrivial behavior is driven by the long-range physics 
associated with the pion exchange and is, in principle, testable in lattice QCD. 

While our analysis establishes correlations between NN low-energy parameters, it does not provide extrapolation 
of lattice-QCD results to the physical point. The possibility of carrying out chiral extrapolations of NN scattering 
observables and the deuteron binding energy by extending the chiral EFT analysis of Ref. m beyond the LO will be 
discussed in a subsequent publication [60] . 








16 


Acknowledgments 


We are grateful to Silas Beane, Dean Lee, Ulf-G. Meifiner and Martin Savage for useful discussions on the topics 
related to this work. We are also grateful to Ulf-G. Meifiner for a careful reading of the manuscript and valuable 
suggestions. One of the authors (EE) thanks the organizers of the YITP-T-14-03 workshop on “Hadrons and Hadron 
Interactions in QGD” at Yukawa Institute for Theoretical Physics, Kyoto University and the organizers of the Bethe 
Eorum on Methods in Lattice Eield Theory, Bethe Center for Theoretical Physics, Bonn University, where a significant 
part of this work was carried out. This work was supported by the European Community-Research Infrastructure 
Integrating Activity “Study of Strongly Interacting Matter” (acronym HadronPhysics3, Grant Agreement n. 283286) 
under the Seventh Eramework Programme of EU, the ERG project 259218 NUCLEAREFT, and by the Georgian 
Shota Rustaveli National Science Foundation (grant 11/31). 


[1] A. Bulgac, G. A. Miller and M. Strikman, Phys. Rev. C 56, 3307 (1997) nucl-th/9708045 . 

[2] E. Braaten and H. W. Hammer, Phys. Rev. Lett. 91, 102002 (2003) nucl-th/0303038 . 

[3] E. Epelbanm, H.-W. Hammer, U.-G. Meifiner and A. Nogga, Enr. Phys. J. C 48, 169 (2006) hep-ph/0602225 . 

[4] P. F. Bedaque, T. Lun and L. Platter, Phys. Rev. C 83, 045803 (2011) arXiv:1012.3840 [nucl-th]]. 

[5] J. G. Berengnt, E. Epelbanm, V. V. Flambanm, G. Hanhart, U.-G. Meifiner, J. Nebreda and J. R. Pelaez, Phys. Rev. D 
87, no. 8, 085018 (2013) |arXiv:1301.1738 [nucl-th]]. 

[6] U.-G. Meifiner, Sci. Bull. (2015) 60(l):43-54 arXiv: 1409.2959 [hep-th]]. 

[7] F. Hoyle, Astrophys. J. Suppl. 1, 121 (1954). 

[8] H. Oberhummer, A. Csoto and H. Schlattl, Nucl. Phys. A 689, 269 (2001) nucl-th/0009046 . 

[9] H. Schlattl, A. Heger, H. Oberhummer, T. Rauscher and A. Gsoto, Astrophys. Space Sci. 291 , 27 (2004) astro-ph/0307528 . 

[10] E. Epelbanm, H. Krebs, T. A. Lahde, D. Lee and U.-G. Meifiner, Phys. Rev. Lett. 110 , no. 11, 112502 (2013) 
arXiv:1212.4181 [nucl-th]]. 

[11] E. Epelbanm, H. Krebs, T. A. Lahde, D. Lee and U. G. Meifiner, Enr. Phys. J. A 49 , 82 (2013) arXiv: 1303.4856 [nucl-th]]. 

[12] D. Lee, Prog. Part. Nucl. Phys. 63, 117 (2009) arXiv:0804.3501 [nucl-th]]. 

[13] E. Epelbanm, H. W. Hammer and U.-G. Meifiner, Rev. Mod. Phys. 81, 1773 (2009) arXiv:0811.1338 [nucl-th]]. 

[14] R. Machleidt and D. R. Entem, Phys. Rept. 503, 1 (2011) arXiv:1105.2919 [nucl-th]]. 

[15] E. Epelbanm, H. Krebs and U.-G. Meifiner, arXiv: 1412.4623 [nucl-th]. 

[16] D. R. Entem, N. Kaiser, R. Machleidt and Y. Nosyk, Phys. Rev. G 91 , no. 1, 014002 (2015) arXiv: 1411.5335 [nucl-th]]. 

[17] E. Epelbanm, H. Krebs and U.-G. Meifiner, arXiv:1412.0142 [nucl-th]. 

[18] S. R. Beane, P. F. Bedaque, M. J. Savage and U. van Kolck, Nucl. Phys. A 700 , 377 (2002) nucl-th/0104030 . 

[19] S. R. Beane and M. J. Savage, Nucl. Phys. A 717, 91 (2003) |nucl-th/0208021 . 

[20] E. Epelbanm, U.-G. Meifiner and W. Glockle, Nucl. Phys. A 714, 535 (2003) nucl-th/0207089 . 

[21] S. R. Beane and M. J. Savage, Nucl. Phys. A 713, 148 (2003) hep-ph/0206113 . 

[22] E. Epelbanm, U.-G. Meifiner and W. Glockle, nucl-th/0208040 

[23] J. Soto and J. Tarrus, Phys. Rev. G 85, 044001 (2012) arXiv:1112.4426 [nucl-th]]. 

[24] N. Barnea, L. Contessi, D. Gazit, F. Pederiva and U. van Kolck, Phys. Rev. Lett. 114 , no. 5, 052501 (2015) arXiv:1311.4966 
[nucl-th]]. 

[25] V. V. Flambaum and R. B. Wiringa, Phys. Rev. G 76, 054002 (2007) arXiv:0709.0077 [nucl-th]]. 

[26] J. Mondejar and J. Soto, Enr. Phys. J. A 32, 77 (2007) nucl-th/0612051 . 

[27] D. B. Kaplan, M. J. Savage and M. B. Wise, Nucl. Phys. B 534, 329 (1998) nucl-th/9802075 . 

[28] T. D. Gohen and J. M. Hansen, Phys. Rev. G 59 , 13 (1999) nucl-th/9808038 ; Phys. Rev. C 59 , 3047 (1999) nucl- 
th/9901065. 

[29] S. Fleming, T. Mehen and 1. W. Stewart, Nucl. Phys. A 677, 313 (2000) nucl-th/9911001 . 

[30] E. Epelbanm, U.-G.-Meifiner, W. Glockle, G. Elster, Phys. Rev. G 65, 044001 (2002). 

[31] E. Epelbanm and J. Gegelia, PoS GD 12, 090 (2013) arXiv:1301.6134 [nucl-th]]. 

[32] E. Epelbanm and J. Gegelia, Phys. Lett. B 716, 338 (2012) arXiv: 1207.2420 [nucl-th]]. 

[33] V. G. Kadyshevsky, Nucl. Phys. B 6, 125 (1968). 

[34] M. Fukugita, Y. Kuramashi, M. Okawa, H. Mino and A. Ukawa, Phys. Rev. D 52, 3003 (1995) hep-lat/9501024 . 

[35] S. R. Beane, P. F. Bedaque, K. Orginos and M. J. Savage, Phys. Rev. Lett. 97 , 012001 (2006) hep-lat/0602010 . 

[36] S. R. Beane et al. [NPLQCD Gollaboration], Phys. Rev. D 85, 054511 (2012) arXiv:1109.2889 [hep-lat]]. 

[37] S. R. Beane et al. [NPLQCD Collaboration], Phys. Rev. D 87, no. 3, 034506 (2013) arXiv:1206.5219 [hep-lat]]. 

[38] T. Yamazaki, K. i. Ishikawa, Y. Kuramashi and A. Ukawa, Phys. Rev. D 86, 074514 (2012) arXiv: 1207.4277 [hep-lat]]. 

[39] T. Yamazaki et al. [PACS-CS Collaboration], Phys. Rev. D 81, 111504 (2010) arXiv:0912.1383 [hep-lat]]. 

[40] S. R. Beane et al. [NPLQCD Collaboration], Phys. Rev. C 88, no. 2, 024003 (2013) |arXiv:1301.5790 [hep-lat]]. 




17 


[41] T. Inoue et al. [HAL QCD Collaboration], Nucl. Phys. A 881, 28 (2012) |arXiv:1112.5926 [hep-lat]]. 

[42] T. Yamazaki, K. i. Ishikawa, Y. Kuramashi and A. Ukawa, arXiv:1502.04182 [hep-lat]. 

[43] S. R. Beane, Nuclear Structure from First Principles, talk at the HHIQCD Workshop, Feb 15 - Mar 21, 2015, Yukawa 
Institute for Theoretical Physics, Kyoto University, Japan 

[44] J. V. Steele and R. J. Furnstahl, Nucl. Phys. A 645, 439 (1999) nucl-th/9808022 . 

[45] E. Epelbaum and J. Gegelia, Eur. Phys. J. A 41, 341 (2009) arXiv:0906.3822 [nucl-th[]. 

[46] E. Epelbaum and J. Gegelia, PoS CD 09, 077 (2009). 

[47] E. Epelbaum, arXiv:1001.3229 [nucl-th]. 

[48] A. M. Gasparyan, M. F. M. Lutz and E. Epelbaum, Eur. Phys. J. A 49, 115 (2013) arXiv:1212.3057 [nucl-th]]. 

[49] J. A. Oiler, arXiv:1402.2449 [nucl-th]. 

[50] J. M. Blatt and J. D. Jackson, Phys. Rev. 76, 18 (1949). 

[51] H. A. Bethe, Phys. Rev. 76, 38 (1949). 

[52] B. Midya, J. Evrard, S. Abramowicz, O. L. R. Suarez and J. M. Sparenberg, arXiv:1501.04011 [quant-ph]. 

[53] M. Pavon Valderrama and E. Ruiz Arriola, Phys. Rev. C 72, 044007 (2005). 

[54] H. van Haeringen and L. P. Kok, Phys. Rev. A26, 1218 (1982). 

[55] A. M. Badalian, L. P. Kok, M. 1. Polikarpov, and Y. A. Simonov, Phys. Rept. 82, 31 (1982). 

[56] E. Epelbaum and U.-G. Meifiner, Few Body Syst. 54, 2175 (2013) nucl-th/0609037 . 

[57] A. Gezerlis, 1. Tews, E. Epelbaum, S. Gandolfi, K. Hebeler, A. Nogga and A. Schwenk, Phys. Rev. Lett. Ill, no. 3, 032501 
(2013) arXiv: 1303.6243 [nucl-th]]. 

[58] A. Gezerlis, 1. Tews, E. Epelbaum, M. Freunek, S. Gandolfi, K. Hebeler, A. Nogga and A. Schwenk, Phys. Rev. C 90, no. 
5, 054323 (2014) arXiv: 1406.0454 [nucl-th]]. 

[59] G. P. Lepage, nucl-th/9706029 

[60] V. Baru, E. Epelbaum, A. A. Filin, and J. Gegelia, “Chiral extrapolations of the NN low-energy parameters,” in preparation. 

[61] V. G. J. Stoks, R. A. M. Klomp, C. P. F. Terheggen and J. J. de Swart, Phys. Rev. C 49, 2950 (1994) |nucl-th/9406039 . 

[62] J. M. Blatt and L. G. Biedenharn, Rev. Mod. Phys. 24, 258 (1952). 

[63] E. Epelbaum, A. M. Gasparyan, J. Gegelia and H. Krebs, arXiv:1501.01191 [nucl-th]. 

[64] J. J. de Swart, C. P. F. Terheggen and V. G. J. Stoks, nucl-th/9509032 

[65] M. Pavon Valderrama and E. Ruiz Arriola, Phys. Rev. C 74, 054001 (2006) [nucl-th/0506047 . 

[66] S. Elhatisari, S. Konig, D. Lee and H.-W. Hammer, Phys. Rev. A 87, no. 5, 052705 (2013) |arXiv:1303.5261 [physics.atom- 
ph]]. 

[67] S. Diirr et al. [Budapest-Marseille-Wuppertal Gollaboration], Phys. Rev. D 90, no. 11, 114504 (2014) arXiv: 1310.3626 
[hep-lat]]. 

[68] G. Alexandrou, M. Gonstantinou, V. Drach, K. Jansen, G. Kallidonis and G. Koutsou, PoS LATTIGE 2013, 292 (2014) 
arXiv:1312.2874 [hep-lat]]. 

[69] C. Alexandrou, M. Gonstantinou, K. Hadjiyiannakou, K. Jansen, C. Kallidonis and G. Koutsou, arXiv:1411.3494 [hep-lat]. 

[70] R. Horsley, Y. Nakamura, A. Nobile, P. E. L. Rakow, G. Schierholz and J. M. Zanotti, Phys. Lett. B 732, 41 (2014) 
arXiv:1302.2233 [hep-lat]]. 



